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1. Introduction 

One common feature in New Physics models is the conservation (or near conservation) 
of a new quantum number. Often it is associated with a parity symmetry, like R-parity 
in super symmetric models or T-parity in Little Higgs ones. Such conserved parity-like 
symmetries serve two basic model-building purposes: firstly, they forbid odd-parity tree 
level corrections to electroweak precision observables, and secondly, they make the lowest 
lying odd-parity state stable, thus providing a possible dark matter candidate. The new 
charge may alternatively come from a continuous symmetry, a global symmetry or a new 
gauge symmetry, for example, and still fulfill the same purposes. 

Regardless of the specific model realization, we can imagine that a new conserved 
quantum number is discovered at LHC. 

In this article, we wish to take some first steps towards addressing a general phe- 
nomenological question: if a new apparently conserved charge should be discovered, is it 
possible to determine experimentally whether it arises from a discrete, a global or a gauge 
symmetry? Specifically, is it possible to determine whether it is the source of a new field? 
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In principle, a continuous symmetry has additive quantum number conservation 
whereas a discrete one has multiphcative conservation. To distinguish between gauge and 
global symmetries one could look for Gauge bosons, for Goldstone bosons and in general 
at the particle spectrum. The new sector however may be "hidden". That is, the carriers 
of the new symmetry, for the two basic reasons mentioned above, could lie entirely within 
the new sector and be neutral under (or have very weak) SM interactions. Indeed, a new 
unbroken symmetry would have to be invisible or else it would have already been found. 
Thus any radiation or other dynamic phenomena associated with it would be invisible to 
SM matter. 

If the charge did radiate in the new sector, would we be still be able to observe indirectly 
the effects of the hidden radiation? How would the kinematic distributions of the visible 
particles be affected? Could we extract information from these kinematic distributions 
about the dynamics within the hidden sector? Could one distinguish Abelian from non- 
Abelian gauge groups, study the different particle (or unparticle) contents or measure the 
strength of the couplings? This could lead to a better understanding of the higher-energy 
dynamics, the ultraviolet completion of the theory, the symmetries involved, and possibly 
even the mechanism by which they are broken. 

The ideal terrain to begin to explore these effects is Hidden Valley models Q. We 
extend this name to the class of models satisfying the following criteria. First, there must 
be a new light hidden sector (the valley), decoupled from the visible SM one, that has 
not yet been discovered because of some barrier. This can be an energy barrier or of 
another nature, e.g. symmetry-forbidding tree-level couplings. Second, the decoupling of 
the new hidden valley sector from the visible SM one must happen at relatively low energies, 
around the TeV scale, in such a way that the cross sections for Standard visible particles 
disappearing into the hidden sector (and vice versa) are small enough to evade the current 
experimental limits, and yet large enough to be observable at LHC. These experimental 
limits are of course model dependent, as we will discuss in sections ^ and ^ 

Typically, the valley particles " v-particles" are charged under a valley group and 
neutral under the SM group Gsm-, and the SM particles are neutral under G„. In order to 
have interactions between the two sectors there has to be a "communicator" which couples 
to both SM and valley particles. A common choice is to have a coupling via a Z' or via 
loops of heavy particles carrying both GgM and Gy charges. 

Examples of Hidden Valleys can be found in many models, such as String Theory Q, 
Twin Higgs models [|, folded SUSY |, |6|, and Unparticle models H §• 

Hidden Valley scenarios can naturally provide candidates for Dark Matter and can 
easily fit cosmological constraints. Just to give an example, in Q the v-interactions ensure 
that all particles efficiently decay to the lightest mesons. These mesons are allowed to 
annihilate to neutral vrJJs, which can then tunnel back into the SM. So long as the lifetime 
of the vrj^ is r <C 1 sec, the number of vrs left will decay exponentially before big-bang 
nucleosynthesis . 

The reason why these scenarios are ideal to study the effects of radiation is the large 
disparity in the masses of the communicators and the v-particles. Typically, the commu- 
nicator has a mass around the decoupling scale, say the TeV scale, while the v-particle 
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mass may be as low as 1-10 GeV. If both the communicator and the v-particle are charged 
under a new gauge they wih radiate gauge bosons, and the larger their mass ratio the 
larger the amount of phase space available for the radiation, both in the normal and in the 
hidden sector. Thus if any effect at all are to be observable, it would be in this kind of 
scenario. 

We have devised a Hidden Valley toy model to tackle the issue, and have implemented 
it in the Pythia 8 Monte Carlo event generator Q . The implementation allows for different 
valley flavour contents, particle masses, gauge groups and v- gauge couplings. In this way 
one may accommodate a range of different Hidden Valley scenarios. 

MC event generators offer flexible approaches to model radiation and parton shower 
evolution in great detail. One new central feature in the Pythia 8 implementation is 
the "competition" mechanism between the hidden and the SM radiation, which is imple- 
mented as an "interleaved" shower, wherein different kinds of emissions, SM and hidden, 
can alternate if viewed in terms of a common shower evolution scale. As a consequence, 
subsequent emissions in the visible sector, of gluons or photons, will then tend to have a 
lower energy than they would have had, had the hidden radiation not been there. This is 
the key mechanism whereby we gain access to the information about the radiation in the 
hidden sector. 

The intention of this article is not a full-fledged experimental analysis of how a new 
sector should be discovered and explored, neither with respect to potential background pro- 
cesses nor to detector-specific capabilities — since our implementation is publicly available, 
we safely leave it to the experimental community to assess. What we want to ascertain 
here is if there are observable signals of hidden valley radiation at all, at the simple parton 
and hadron levels. 

It is not trivial to decide which visible particle kinematic distributions one should 
study to reveal valley radiation effects and to discriminate between different models. For 
instance, at an e^e~ collider the rise of the communicator pair-production cross section 
near threshold could allow to determine its spin, and thereafter the absolute size of the cross 
section could suggest the presence of new "colour" factors — recall that the pair-production 
of particles in the fundamental representation of a new SU{N) group gives a factor in 
the cross section. Such measurements would not directly probe the hidden sector, however: 
they would not reveal whether a new group is gauged, or what is the coupling strength in 
it. For a hadron collider, like the LHC, the uncertainty in the event-by-event subcollision 
energy \/I undermines analyses solely based upon the value of cross-section. The best 
strategy is thus to complement cross-section with invariant mass measurements and the 
study of other boost-invariant quantities (for a recent review see the proceedings [^]). 

This is the reason why we choose to study MT2 [15| distributions, which give relations 
between communicator and v-particle mass. These observables are specifically designed 
to be boost-invariant and to deal with BSM models in which more than one particle es- 
capes detection, such as in our toy model. But we also study "hidden observables", like 
the invariant mass distribution of a hidden particle together with its associated hidden 
radiation. 

The effects of the hidden radiation on these distributions and how much one may 



-3- 



observe depends on the details of the scenario considered, of course, but also depend 
heavily upon the collider type considered, on its center of mass energy, and on its integrated 
luminosity L. We consider two different LHC scenarios, one for the early data (the first 
18 to 24 months at 7 TeV with an expected integrated luminosity £ = 1 fb^^) and one for 
later data (-y/i = 14 TeV and integrated luminosity C = 100 fb^^). The conclusions in the 
two cases will be quite different. For e^e~ collisions we will mainly refer to an ILC at 800 
GeV, though we mention CLIC production cross sections at 3 TeV. 

The paper is organized as follows: we give a first overview of the various valley scenarios 
in the literature, in section |2|, and then describe the tools that are now available in Pythia 
in section ^ In section ^ we explain the model and the main features of interleaved SM and 
hidden radiation. Finally, in the last two sections we study the effects of hidden radiation 
on collider phenomenology. In section ^ we discuss the e^e~ case and in section |6| the LHC 
one. 

2. Hidden Valley scenarios 

As mentioned in the introduction a Hidden Valley is a light hidden sector, consisting of 
particles which, depending on the model, might have masses as low as 10 GeV. The detailed 
spectrum of the v-particles and their dynamics within the hidden valley depends upon the 
valley gauge group G^,, the spin and number of particles present in the theory, and the 
representation they belong to. 

The effects of the hidden sector on the visible particle spectra will depend upon the way 
the hidden sector communicates with the SM, whether it is via a Higgs, multiple Higgses, 
a Z\ heavy sterile neutrinos or via loop of heavy particles charged under both SM and 
valley gauge interactions. 

We would like to give a panoramic view of the different Hidden Valley scenarios without 
going into details and to underline those features that may be simulated with the new tools. 

The simplest possibility is a QCD-like scenario, with a strong coupling constant, which 
may run like the QCD coupling does, with QCD-like hadronization generating valley pions, 
v-r/s, v-ifs, v-nucleons etc. The Standard Model 5C/(3)c x SU{2) x U{1) sector could 
couple ultra-weakly with the hidden SU{N) sector via a neutral Z'. This scenario was 
investigated by Strassler and Zurek with tools analogous to the ones used to simulate 
QCD H]. It displays some rather startling features. For instance, a v-vr could have a 
displaced decay in the muon spectrometer in the ATLAS detector, resulting in a large 
number of charged hadrons traversing the spectrometer, or it could decay in the hadronic 
calorimeter producing a jet with no energy deposited in the electromagnetic calorimeter 
and no associated tracks in the inner detector. Experimental studies for these scenarios 
are currently under way, by the DO, CDF, LHCb, ATLAS and CMS collaborations. 

Typical hidden valley-like signatures appear also in Unparticle models with mass gaps 
0. These models display a conformal dynamic above the mass gap, and a hidden valley 
behaviour when the conformal symmetry is broken. Regardless of the dynamics above 
the mass gap, whether it is strongly coupled or weakly coupled, the signatures are similar 
to the ones mentioned in the previous scenario (displaced vertices and missing energy 
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signals). This is because only the lower energy states, light stable hidden hadrons can 
decay back into Standard Model particles. The higher energy states, be they narrow 
resonances or a continuum of resonances, decay rapidly to these lower light stable hadrons. 
As for the previous scenario, a parton shower is a key tool to study these models, not so 
much to determine the phenomenology qualitatively, but because it the only element of the 
hidden dynamics which is sensitive to the higher-energy conformal (or next-to-conformal) 
dynamics. The conformal dynamics will be reflected in the parton shower evolution, which 
can be rather different from the regular QCD one, especially in theories with a strong 
dynamics above the mass gap. 

There are of course many other Hidden- Valley related models, such as Quirky models 
Q, just to give an example, in which the parton shower evolution does not play the key 
role it does in the previous cases. Typically their phenomenology is better captured in 
terms of string dynamics and string fragmentation. 

In this paper we do not address the issue of string fragmentation or hadronization. 
Our main focus is on the parton shower, as this best captures the nature of the hidden 
radiation. 

The model we built to investigate the existence of a this new radiation exploits but 
a few features common in many hidden valley scenarios: the presence of a new unbroken 
hidden gauge group, of a heavy communicator, charged under both SM and hidden sector 
gauge group and decaying into a visible and an invisible light particle, charged only under 
the new gauge group. These characteristics fit many Hidden Valley models, we however 
make an additional assumption, which is that the production cross sections^ should be 
large enough for the effects of the hidden radiation to be discernable. This model was then 
implemented in the Pythia event generator. Notice however, that the shower mechanism 
we implemented is rather different from the ones mentioned above, as we will discuss in 
the next section. 

3. Monte Carlo Tools in PYTHIA 8 

In order to allow detailed studies of a set of scenarios, the models have been implemented 
in the Pythia event generator, and will be publicly available from version 8.140 onwards. 

3.1 Particle content 

For simplicity we assume that the HV contains either an Abelian U{1) or a non-Abelian 
SU{Nc) gauge group, with spin 1 gauge bosons. The former group could be unbroken or 
broken, while the latter always is assumed unbroken. Casimir constants could be general- 
ized to encompass other gauge groups, should the need arise, but for now we do not see 
that need. The gauge bosons are called 7^ and g^j, respectively. 

A particle content has been introduced to mirror the Standard Model flavour structure. 
These particles, collectively called F^j, are charged under both the SM and the HV symmetry 
groups. Each new particle couples flavour-diagonally to a corresponding SM state, and has 

""^ We will discuss the production cross sections in section W. 
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Table 1: The allowed particle content in the HV scenarios, with their SM partners, where relevant. 
The code is an integer identifier, in the spirit of the PDG codes, but is not part of the current 
Ainsler:2008zzb standard. 

the same SM charge and colour, but in addition is in the fundamental representation of 
the HV colour, see Table |. Their masses and widths can be set individually. It would also 
be possible to expand the decay tables to allow for flavour mixing. 

These particles can decay to the corresponding SM particle, plus an invisible, massive 
HV particle qy, that then also has to be in the fundamental representation of the HV colour: 
— )• fq^. The notation is intended to make contact with SM equivalents, but obviously 
it cannot be pushed too far. For instance, not both Fy and q^ can be fermions. We allow 
the Fy to have either of spin 0, 1/2 and 1. Currently the choice of q^ spin is not important 
but, for the record, it is assumed to be spin 1/2 if the F^ is a boson and either of spin 
and 1 if F^ is a fermion. 

3.2 Production processes 

The HV particles have to be pair-produced. The production processes we have implemented 
are the QCD ones, gg — )• QvQv and qq — )• QvQv^ fo^' the coloured subset Qv of F^ states, 
and the electroweak // — )• ^* jZ^ — )• for all states. All of them would contribute at a 

hadron collider, but for a lepton one only the latter would be relevant. Each process can 
be switched on individually, e.g. if one would like to simulate a scenario with only the first 
Fy generation. 

Note that pair production cross sections contain a factor of N^, with A'^r; = 1 for an 
C/(l) group, for the pair production of new particles in the fundamental representation of 
the HV gauge group, in addition to the ordinary colour factor for Q^. Other things equal, 
this could be used to determine A^c from data, as already discussed. For the case of a spin 
1 Fy it is possible to include an anomalous magnetic dipole moment, k^\. 

The spin structure of the F^ — )• /g^, decay is currently not specified, so the decay is 
isotropic. Also the Yukawa couplings in decays are not set as such, but are implicit in the 
choice of widths for the F^ states. 

The kinematics of the decay is strongly influenced by the q^ mass. This mass is almost 
unconstrained, and can therefore range from close to zero to close to the F^ masses. We 
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will assume it is not heavier than them, however, so that we do not have to consider the 
phenomenology of stable Fy particles. 

3.3 Parton showers 

Both the and the g„ can radiate, owing to their charge under the new gauge group, i.e. 
Fy y Fqj^qj and qqj )■ Qv^v fo^ ^ ^(-'-) S^^o^P) s-nd F.(, y F^jg^j and Qvgv for a SU{Nc) 

one. In the latter case also non-Abelian branchings — )■ gvOv are allowed. Currently both 
and gv are assumed massless, but a broken U{1) with a massive is foreseen. 
These showers form an integrated part of the standard final-state showering machin- 
ery. Specifically, HV radiation is interleaved with SM radiation in a common sequence of 
decreasing p±. That is, at the stage before the F„'s decay, they may radiate g, 7 and 'yv/gv, 
in any order. For the z'th emission, the p±_ evolution starts from the maximum scale given 
by the previous emission. The overall starting scale p_Lo is set by the scale of the hard 
process. Thus the probability to pick a given p_\_ takes the form 

dp ^ AIPqcd ^ dpQED ^ dPiivA / /"f^-i /dpQCD ^ dT^QED ^ dpHv \ \ 

(3.1) 

where the exponential corresponds to the Sudakov form factor. Implicitly one must also 
sum over all partons that can radiate. 

To be more precise, radiation is based on a dipole picture, where it is a pair of partons 
that collectively radiates a new parton. The dipole assignment is worked out in the limit 
of infinitely many (HV or ordinary) colours, so that only planar colour flows need be 
considered. Technically the total radiation of the dipole is split into two ends, where one 
end acts as radiator and the other as recoiler. The recoiler ensures that total energy and 
momentum is conserved during the emission, with partons on the mass shell before and 
after the emission. In general the dipoles will be different for QCD, QED and HV. 

To take an example, consider qq — > QvQv, which proceeds via an intermediate s- 
channel gluon. Since this gluon carries no QED or HV charge it follows that the QvQv pair 
forms a dipole with respect to these two emission kinds. The gluon does carry QCD octet 
charge, however, so QvQv do not form a QCD dipole. Instead each of them is attached to 
another parton, either the beam remnant that carries the corresponding anticolour or some 
other parton emitted as part of the initial-state shower. This means that QCD radiation 
can change the invariant mass of the QvQv system, while QED and HV radiation could 
not. When a 7 or 7^, is emitted the dipole assignments are not modified, since these bosons 
do not carry away any charge. A g oi gv would, and so a new dipole would be formed. For 
QCD the dipole between Qv and one beam remnant, say, would be split into one between 
the Qv and the g, and one further from the g to the remnant. For HV the QvQv dipole 
would be spit into two, Qvgv and gvQv As the shower evolves, the three different kinds of 
dipoles will diverge further. 

Note that, in the full event-generation machinery, the final-state radiation considered 
here is also interleaved in p± with the initial-state showers and with multiple parton-parton 
interactions. 
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There is made a clean separation between radiation in the production stage of the 
t, pair and in their respective decay. Strictly speaking this would only be valid when 
the width is small, but that is the case that interests us here. In the decay F^ — >• fq^ the 
QCD and QED charges go with the / and the HV one with g^. For all three interactions 
the dipole is formed between the / and the qv, so that radiation preserves the F^ system 
mass, but in each case only the relevant dipole end is allowed to radiate the kind of gauge 
bosons that goes with its charge. (Strictly speaking dipoles are stretched between the / 
or q^j and the "hole" left behind by the decaying F^. The situation is closely analogous to 
t — )- bW~^ decays.) 

The HV shower only contains two parameters. The main one is the coupling strength 
a„, i.e. the equivalent of a^. This coupling is taken to be a constant, i.e. no running is 
included. 

From a practical point of view it is doubtful that such a running could be pinned down 
anyway, and from a theory point of view it means we do not have to specify the full flavour 
structure of the hidden sector. The second parameter is the lower cutoff scale for shower 
evolution, by default chosen the same dbS for the QCD shower, pj_min — 

0.4 GeV. 

The HV showers are not matched onto higher-order matrix elements for the emissions 
of hard ^v/dv in the production process, and so contain an element of uncertainty in that 
region. For the decay process the matching to first-order matrix elements have been worked 
out for all the colour and spin combinations that occurs in the MSSM 11], and is recycled 
for the HV scenarios, with spin 1 replaced by for non-existing (in MSSM) combinations. 
This means that the full phase space is filled with (approximately) the correct rate. Some 
further approximations exists, e.g. in the handling of mass effects in the soft region. The 
chosen behaviour has been influenced by our experience with QCD, however, and so should 
provide a good first estimate. More than that we do not aim for in this study. 

4. The model: SM and SU{3)^ radiation 

To be specific, in the following we explore two similar Hidden Valley experimental scenarios. 
In the first, the communicator Ey is a spin 1/2 particle charged under both the SM SU{2) x 
U (1) and the valley gauge group SU{3)y. We assume it has the same SM charges an electron 
would have, so it may be pair-produced in e~^e~ collisions, via Z/j*. Under the unbroken 
SU{3)y, it transforms like a 3, so it radiates both 7s and massless hidden valley gluons 
QyS. After the parton shower, the E^ eventually decays into a visible SM electron e and 
an invisible spin valley "quark" q^. This q^ belongs to the fundamental representation of 
SU{3)j} and is not charged under the SM gauge group, so it only radiates g^s. See Fig. ||. 

The key feature is interleaved radiation, already introduced above. In the current 
context it works as follows. Once the Ey has been produced it may radiate a SM 7, say. 
This radiation will subtract energy from the Ey and the following emission, be it another 
SM photon or a valley gluon, will have less phase space to radiate into. In an analogous 
way, assuming a valley gy is emitted next, it subtracts energy from Ey and affects the 
following emissions which, again, could be either visible or invisible. 
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Figure 1: An EyEy pair is produced via Z/j*. Since Ey is charged under both SU{2) x U{1) and 
SU{5)y, it radiates both 7s and gyS. It eventually decays into e and qy. These then each radiate 
into their respective sector. Notice that qy here refers to a spin particle. 




Figure 2: Pair production of hidden valley DyS. Each Dy can radiate ga and gyS, and eventually 
decays Dy — > qyd. The visible d then can radiate further gs and the invisible qy further gyS. 

In the second scenario the communicator between SM sector and Hidden Valley sector 
is a quark-like, spin 1/2 object Q^, belonging to the (3, 3) representation of the gauge group 

vS are pair produced (mostly) via strong interactions (gluon-gluon 
or qq fusion). We choose the scenario in which only one vector-like is produced, the D^. 
This emits massless valley gluons (since the SU{3)v is assumed to be unbroken) and 
these may in turn radiate more g^s. During the shower evolution, both types of gluons are 
radiated until finally each decays into a visible SM d quark an invisible spin valley 

The decays are flavour diagonal, ^ d + q^. 

The SM quark d transforms as a (3, 1) under 5C/(3)c x SU{?>)y, so it radiates only SM 
gluons, while the valley q^ belongs to the (1,3) representation of SU{2>)c x 5C/(3)^, so not 
having any SM color charge, it radiates only g.uS., see figure ^. 

In both scenarios there are just three parameters left to vary: the size of the valley 
coupling constant q„, the masses of the communicator particles Me^ or M/)^ and the mass 
of the valley scalar Mq^ . 

Below, in Table ^, we list the total production cross sections at different colliders: 
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ILC 
(800 GeV) 


CLIC 
(3 TeV) 




LHC 

(7 TeV) 


LHC 
(14 TeV) 


Me, = 300 GeV 


398 fb 


44 fb 


Md, = 300 GeV 


1.39 • 10^ fb 


1.04 • 10^ fb 


Me, = 500 GeV 




41 fb 


Md, = 500 GeV 


654 fb 


7.27 • 10^ fb 


Me, = 1 TeV 




32 fb 


Md, = 1 TeV 


3.21 fb 


124 fb 



Table 2: The order of magnitude of the total production cross sections, in /6, at ILC (via Z/7*), 
LHC (via qq or gg fusion) with ^/s = 7 TeV and 14 TeV, for various values of the communicator 
mass. The spin of the communicator is assumed to be 1/2. 
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Figure 3: The spin dependence of the EyEy production cross section at e+e^ colliders for Ey spin 
1/2, and 1. Me, = 250 GeV, the mass spread Te, = 2 GeV, Mg„ = 50 GeV. The spin curve 
has been scaled by a factor 30. 

e+e- with = 800 GeV or Vs = 3 TeV, and LHC with ^/s = 7 TeV or Vs = 14 TeV for 
some typical Me, , Md, mass values. 

We also show the spin dependence of the EyE^ production cross section at e"'"e~ 
colliders for the three cases: Fy spin and spin 1/2, Fy spin 1/2 and Qy spin or 1, and 
Fy spin 1 and qy spin 1/2, Fig. ^ 

The higher the spin, the larger the cross section. Indeed, the curve corresponding to 
Ey spin has been scaled by a factor 30 to emphasize the similarity in shape with the 
spin 1 curve. Note that the processes proceed through the s-channel exchange of a spin 1 
7*/Z*. Thus the production of a spin 1/2 pair has only a threshold factor /3 from phase 
space, where /3 is the velocity of the produced pair in the j* /Z* decay vertex, while the 
other two have an (approximate) additional factor from helicity considerations. The 
results have again been obtained with a SU{3)y group, and are directly proportional to 
the Nc chosen. Since we would not expect gauge groups with Nc above (some multiple 
of) 30, the conclusion would be that a threshold scan of the cross section could be used 
to determine both the Ey spin and the number of hidden colours, as well as the Ey mass, 
of course. A caveat would be that we have here only considered the Ey gauge production 
mechanism, not the possibility of a significant t— channel Yukawa contribution. The reason 
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Figure 4: Left: the effective invariant mass distribution for the at LHC for = 0.1,0.4,0.8. 
Afc,, = 300 GeV, nominal g„ mass Mq^ — 10 GeV and y/s ~ 7 TeV. Right: the effective invariant 
mass distribution for the qi, at LHC for — 0.1,0.4,0.8. A/^i^, = 1 TeV, nominal qy mass 



Mg^ = 10 GeV and 
bare mass rUo 



14 TeV. Notice how the mean value of the distribution shifts from the 
TTiD^ GeV as the coupling constant grows. 



10 GeV towards m^^'^'^ 



for not including this production channel is that it would imply a large decay width for 
the , which would give additional large and model dependent effects to the cross section 
around threshold (see below). 

The experimental constraints on these two types of setup are similar to the ones for 
New Charged Leptons and Leptoquark production^. For the New Charged Leptons the 



PDG |1C] gives the lower bound mi± > 100.8 GeV. For scalar and vector Leptoquark states 
we use the direct limits coming from leptoquark pair production and subsequent LQ — )• z^g 
decay searches |11]. For the all generation search rriLQ > 136 GeV in the scalar case and 



rriLQ > 200 GeV in the vector case, where one assumes the branching ratio B{vq) = 1. 
Bounds on the leptoquark mass for the third generation decaying into vb are more stringent 
0, m^Q > 229 GeV. 

We view these boundaries as simply indicative of the mass range we should contemplate 
and chose masses that are well beyond these boundaries. Our studies are in any case not 
critically dependent on them. 

We assumed the communicators to be massive, Me^ in the range [250,300] GeV for 
the ILC case, Md, in the range [300, 500] GeV for the LHC 7 TeV run, and [0.5,1] TeV 
for the 14 TeV run. The hidden scalar is taken to be light, Mq^ as light as 10 GeV. The a„ 
parameter is allowed to vary over a wide range, and results shown for interesting values. 

Depending on the size of these parameters, the effects of the radiation on the lepton or 
quark kinematic distributions can be significant, as we will show in the next two sections. 
But whether these effects will be observable strongly depends on the statistics at hand. 
We assume an integrated luminosity L = 200 fh~^ for the ILC, 1 fb^^ for the 7 TeV LHC 
run and 100 fb"^ for the 14 TeV one. 

The hidden radiation affects the visible particle kinematics through two mechanisms. 
The first one, interleaved radiation from the or Dy, we already discussed. The second is 



^We do not make use of the limits coming from the DO or CDF q ^ j + searches, because these 
depend upon the chosen mSUGRA scenario. 
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radiation off the qy after E^/Dy decay. Tliis causes the invariant mass for the system made 
out of Qy and the radiated gyS to be larger than the on-shell qy mass. This invariant mass can 
be viewed as the off-shell mass with which the Qy is produced at the decay vertex. In Fig. ^ 
one may see the effects of the valley radiation on the off-shell mass, for the Dy — y d -\- Qy 
case. As the hidden valley coupling Oy increases, the Qy radiates more and more into the 
hidden sector, and the mean value of the distribution shifts from Mg^ towards M^)^. At 
the same time more and more energy is subtracted from the visible recoiling particle, in 
this case the d and its system of emitted gluons. 

The effect is more obvious when the mass difference Mj)^ — M^^ is large, since more 
phase space is available for the radiation. 

The size of the deviations induced by these two combined mechanisms is very much 
dependent on the collider, as we already stressed above. In the next two sections we will 
discuss the various cases separately. 

5. Effects of SU{3)c radiation at e^e colliders 

We begin by studying the e~^e~ — )• 'y* /Z — )• EyEy, scenario, which allows for many simpli- 
fications compared to the quark case, and therefore offers a convenient warmup. For the 
ILC with ^/s = 800 GeV and an assumed integrated luminosity of L = 200 fb~^ per year 
an Me^ = 300 GeV translates into about 80000 EyEy pairs. 

5.1 Collisions in the center-of-mass frame 

To illustrate the principles, as a very first step we will neglect bremsstrahlung and beam- 
strahlung. We then only need to consider two types of interactions, electromagnetic and 
valley SU{3)y radiation in the final state, with coupling constants a and a„. No fragmen- 
tation or hadronization need to be taken into account. 

Since the center of mass (CM) of the collision is at rest, there is a clean relationship 
between the mass of the hidden valley qy, Mg^, and that of the communicator Me^- In the 
absence of radiation (hidden or standard), this can be inferred from the distribution of the 
energy of the emitted electrons, in particular from the upper endpoint of this distribution, 
describing the electron maximum energy. This is obtained when the electron is emitted in 
the same direction as the Ey is moving in, with the Qy in the opposite direction. One may 
use this maximization condition to derive the relationship between Mg^ and Me^- 

In the rest frame of the Ey, neglecting the electron mass, 

Pe. = (Me., 0,0,0) , 

( Ml + M2 Ml - M2 \ 

^^"-(^ 2Me. 2Ml ) 

Ml -M2 Ml -M^ 

^ 2Me. ' ' ' 2Me. ' 
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(5.1) 



Assuming the boost to the CM rest frame is at an angle 9 with respect to the e direction 
in the rest frame, the electron energy will be given by 




= 7(Ee + /3|Pe| cos 0) = ^ ( 1-^1 I 1 + ^1-1^ COS0 I , (5.2) 

where cos9 = ±1 gives the upper and lower edge of the energy spectrum. If the decay is 
assumed isotropic, dP/dcos0 = constant, the electron energy spectrum is flat between the 
limits. 

So if one can measure the maximum and minimum energy E'^^ one may solve for Me^ 
and Mq^. Fig. ^ shows the energy distribution E^ of the electrons produced with and 
without hidden radiation. In the latter case the spectrum is shifted to lower values, as 
the hidden sector takes a bigger fraction of the available energy, by radiation off both the 
E^ and the q^. The endpoints remain the same, as there is always a fraction of events 
where radiation is negligible. As we have assumed a modest width of 1 GeV for the Ey 
there is a tiny tail beyond the expected edge. (We could cope with a wide range of widths, 
but have picked values in the GeV range, so that the possibility of a Breit-Wigner-shaped 
mass broadening is not overlooked, while still maintaining a credible simulation in terms of 
resonance diagrams only.) The key point to observe, however, is how the upper "shoulder" 
is softened by the hidden radiation. Thereby a precision measurement of this region would 
offer a direct check on the amount of hidden radiation. At the lower end, QED cascades 
such as e~ — )• e~7 — )• e~e'^e~ contribute to the spectrum, but are easily eliminated if only 
the highest-energy lepton is considered, right side of Fig. ^, or at least only the highest 
two. We should clarify that the electron energy studied in this section includes photons 
emitted near the electron direction, since we here include a Durham "jet" algorithm that 
clusters photons within a 3 GeV pTSva.9/2 distance of the electron. 

Whether and how well one would actually be able to observe these endpoints will be 
very model and detector dependant. Regardless of the background or detector sensitivity, 
we expect that the endpoints of the distribution will have low statistics, given that they 
correspond to extreme kinematical configurations. Most likely, one will need to rely on data 
points in the shoulder region to fit the curve and extrapolate the endpoint E'^ These 
shoulder data points would be the ones most affected by the radiation, so the mass Me^ and 
Mq^ inferred from them would be significantly different when hidden radiation is included. 
On the one hand, the curve corresponding to having valley radiation is always softer than 
the one without, so some mean of the threshold region will give too low an endpoint. On 
the other hand, if one only tries e.g. a linear fit, the shape of the fall-off in the threshold 
region would suggest too high an endpoint. A readiness to include a parametric shape for 
the endpoint region, that takes into account a tuneable radiation contribution, will help 
ensure a better extraction of the relevant mass parameters. 

Notice that the curve corresponding to having valley radiation always lies below the 
one without, implying that the value of -Be,max with the radiation would always be higher 
than the one without. This is a reflection of the two mechanisms we mentioned in the 
previous section: first, the the valley gluons subtract energy from the E^, ultimately 
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Figure 5: the energy distribution of the vis- 
ible electrons. Looking at the upper shoul- 
der of the distribution, the upper (black) 
and lower (red) curves give the prediction for 
Ey — >■ e~qy when hidden valley radiation is 
not or is take into account; in both cases elec- 
tromagnetic radiation is included. Center of 
mass energy Vs = 800 GeV, Me^ = 300 
GeV, = 50 GeV and a„ = 0.05. Number 
of events per 6 GeV bins, luminosity L — 200 
fb^^. The error is purely statistical. 
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Figure 6: The energy distribution of the 
most energetic electron in each event, under 
the same conditions. 



subtracting it from the es, and second, when they are emitted by the q^, they change its 
effective mass, M^^ > Mq^, as one may see in Fig. ^ again subtracting energy from the 
decay e. 

Would it be possible to describe the curves with hidden radiation using a model without 
it, but with different mass parameters Me^ and Mq^7 Fig. |^ shows the effects of changing 
the invisible particle mass Mq^ in model with and without radiation. The "fingerprint" 
of the v-radiation is clear: a softening in the shoulder of the distribution which leaves 
the endpoints fixed. A simple change in the mass parameters of the model without hidden 
radiation (in this case Mq^) changes the endpoints and leaves the sharp drop of the shoulder 
unchanged. 

Notice how so long as the mass difference Me^ — Mq^ > 40 GeV one may always dis- 
tinguish between any two curves with and without radiation. This of course dependent. 

In Fig. ^ instead one can see the dependence of the energy distribution for the 
visible electrons in each event. Notice how even for a coupling as low as = 0.05 the 
effects of the v-radiation on the shoulder region are already sizeable. 

There are some parameter regions (e.g. when the Ey-to-q^ mass splitting is small) where 
the shape of the distribution of the hardest leptons is no longer conclusive in distinguishing 
between a model with and one without valley radiation. In this case one may consider 
other observables which have an "orthogonal" dependence on the valley parameters. We 
studied 77, linearized sphericity S and the number of emitted leptons, for example. There 
is no unique strategy in this case, one must perform a case by case study of the different 
observables in the different parameter regions to determine which one displays the largest 



- 14 - 



0.16 




[GeV] 

Figure 7: We compare the energy distributions of ttie most energetic electron emitted in each 
event, in the case ^/s = 800 GeV. The communicator mass is fixed at Me^ = 200 GeV, while Mq^ 
is allowed to vary between 20 and 190 GcV. The valley gauge coupling is fixed at =0.1 in order 
to isolate the mass dependence. Each endpoint corresponds to two curves, the lower one being the 
one with and the top one the one without hidden valley radiation. 
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Figure 8: av dependence of the energy of the visible electrons. The (black) squared-off curve 
corresponds to the model with no valley radiation, the uppermost if one looks at the shoulder. 
Below it are the curves corresponding to valley coupling constant ay = 0.05 and av = 0.1. Number 
of electrons per 6 GeV bin, = 800 GeV, Me, = 300 GeV and Af,„ = 50 GeV. 



separation between the model with and the model without radiation. As a general rule, 
we found that three observables were normally sufficient to distinguish the two. 
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5.2 Collisions not in the center- of- mass frame: MT2 

We now consider the effect of initial-state radiation (ISR). This causes unobservable radi- 
ation, mainly along the beamline, but also some transverse kicks. Beamstrahlung is highly 
machine-dependent and thus not included, but is purely longitudinal. The methods we will 
introduce to handle bremsstrahlung also automatically handle beamstrahlung with little or 
no degradation of performance, so from now on we will not address the latter specifically. 

For our theoretical studies, in order to avoid the clustering of ISR 7 radiation with the 
leptons coming from the hard interaction, we apply a cut on the r] > 5. The symmetry of 
the system now being cylindrical, we also changed the clustering algorithm to the cylindrical 



fast jet [14 



The major consequence of ISR is that the collision now no longer happens in the CM 
rest frame, with the information connected to the pz, the momentum along the beampipe, 
no longer available. In this case it is convenient to introduce a new variable called Cam- 



bridge MT2, see ^ 



The MT2 variable was invented precisely to treat events in which the new particles are 
pair-produced and then each decay into one particle that is directly observable and another 
particle whose existence may only be inferred from from missing transverse momenta. 

This observable is somewhat inspired by the transverse mass mj- used at hadron col- 
liders to measure the mass of the W boson in the decay W — )• ez^. The neutrino escapes 
detection, its only trace in the detector being missing momentum. In this case one can 
construct the variable 

= 2{E^ - p^- ^^). (5.3) 



Here Et is defined as Et = y + p^, although in this particular case the electron and 
neutrino masses can be neglected, of course. The rn^ variable has the property that 

my < rriyy (5.4) 

If there is enough statistics to ensure that the kinematic configuration corresponding 
to the maximum is hit, this gives a measurement of (a lower bound on) the W mass. 
Analogously, one may build a variable called MT2, with the property that its upper bound 
describes the mass of the communicators, i.e. the particles that were pair-produced. 

Now consider the process described in Fig. |9|. Two particles Y are pair- produced, 
then each of them decays into a visible particle (a or b in the figure) and one that escapes 
detection, called A'^i_2- One cannot use the transverse momentum in this case since there are 
two particles escaping detection, both contributing to the missing transverse momentum 

The observable MT2 is defined as 

MT2 = min [max {my(pra, F^Ti), "^tCptJ), ^T2)}] • (5.5) 

F*Ti+F*T2— Pt 

where ]^ti,T2 are all the possible 2-momenta taken away by the A^s, such that their sum 
gives the observed missing momenta, t^ti+ F^T2 =F^t- 

MT2 coincides with the mass of the communicator Y, i.e. MT2 has a maximum, when 
for both communicator decays the visible and the invisible particles are produced at the 
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UTM 




Figure 9: The MT2 diagram for the short- 
est, simplest decay chain . Two communica- 
tor particles are pair-produced and each de- 
cay into a visible (a or b) and an invisible (A^i 
and N2) particle. Upstream transverse mo- 
mentum (before the decay) is also possible, 
but it must be known. Longitudinal momen- 
tum information is not available. 
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Figure 10: An example of the MT2 method 
applied to the communicator mass Me^ ■ His- 
togram of the MT2 values obtained for a 
communicator mass Me^ = 200 GeV and 
Alq^ and an invisible particle mass M^^, = 
100 GeV. The small tail in the distribution 
is due to the F = 1 GeV spread in the mass 
distribution (more about his in the follow- 
ing). 
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(5.6) 



For the studies in this article we used a particularly simple version of the MT2 algo- 
rithm 1^], the source code of which can be downloaded from 
http:/ /daneel. phyics.uc.davis.edu/ Cheng:2008hk/mt2-l. Ola/test. 



A more sophisticated algorithm is described in 



The simpler method is based on the 



use of "kinematic constraints" |21, 22, 17, 20 1 



P? 



(Pl 



Paf = (P2 +Pb)^ 



My 



Pl +P2 = 
P1+P2 = f 



(5.7) 



In the case of two invisible and at least two visible particles as in Fig. the two 



methods actually coincide [20|. 

The inputs of the MT2 method are rriN, iria, mi,, pf^, p^^ and jS-r- Notice that rua 
and ruf, may change quite substantially from event to event, since they each correspond 
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Figure 11: The effect of valley radiation on the variable MT2 for different a„ — 0.05, 0.1, 0.2 values 
(close-up on the shoulder region) at CLIC. The distributions were obtained assuming a luminosity 
L = 1 fb-\ = 1 TeV and Mq„ = 50 GeV. Left: the effect for a hypothetical Me^ = 200 GeV. 
Right: the effect of valley radiation on MT2 for a hypothetical Me^ — 300 GeV. 



to the invariant masses of the clustered visible particles (in this case the lepton and the 
photons) of each branch. The output of the MT2 method is one single MT2 value per 
event. Fig. IC illustrates a typical use of the MT2 variable. If one histograms the MT2 
values over a large number of events, the upper edge of this distribution gives a lower limit 
on the communicator mass My. 

Whether the event rate in the upper-edge kinematic region defined in eq. ( ^.6D is large 
enough to be able to extract the endpoint MT2™^^ = My, for a given luminosity, depends 
of course on the interactions. Even more than for the energy variable in the previous 
section, it is not unlikely that MT2™'^'^ might have to be extrapolated from points in the 
shoulder region. 

There are many other methods to determine mass relations between the the new par- 
ticles, [17, just to cite some. Some of these are very closely related to MT2, such 
as Ijl^. Some of these require cascade decay chains, or make assumptions about the new 
particles involved in the decay chain being on shell, or require high luminosity. Where this 
information is actually available, one should of course make use of it, |24, 

The assumptions in this study, though, are that each of the identical decay chains 
consists of a single two-body decay and that the integrated luminosity, at least for the 
LHC at 7 TeV study, might be rather low (1 fb~^). These effectively preclude the use of 
many of the above methods. 

In Fig. 11 one may see the effect of the valley radiation on the MT2 distribution for 
different a„ values and communicator mass parameters. The interesting region is again 
represented by the "shoulder" of the distribution. Notice how the amount of invisible 
radiation, and thus the effect on MT2™'^^, increases with Me^ — M^^, analogously to what 
happens in the energy distributions. The size of these effects may be compared with the 



effects coming from ISR, in Fig. 12 



The MT2 distribution might present a tail, due to the Breit-Wigner spread F that 



we allow for, see Fig. 13. As one may see, all the MT2 data points which lie above 
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Figure 13: Left: the (right side of the) Breit-Wigner mass distribution of the communicator E^, 
with an Mq = Me, = 300 GeV and a width = 3 GeV. Right: MT2 distribution for the 
same points. The upper (red) points correspond to considering only EyS having masses within the 
IfnE^ — Me^, \ <^e^ GeV interval. The lower points (green) describe the MT2 points corresponding 
to both E^s having a mass greater than Me^, + ^e,- All other effects have been switched off. 
Number of events per 1 GeV bin, ^/s — 800 GeV, Mq, = 50 GeV, integrated luminosity L — 200 
fb^^. The error is purely statistical. 
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300 GeV correspond to EyS which actually have a larger mass than the 



As we already stated above, niiy, 
the mass of the invisible particle, is 
an input parameter. MT2™'^^ only 
gives one single relation between My 
and TTiTv- Depending upon the de- 
cay chain topology and the presence 
or not of upstream transverse mo- 
mentum (UTM, in the following, may 
come from ISR or from previous de- 
cays), there are different strategies 
to determine both masses simultane- 
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Figure 12: The effects of initial-state radiation on the 
MT2 distribution and on the value of the Me^ mass one 
infers. Number of events per 4 GeV bin, for an e+e^ 
collider with = 1 TeV, L = 1 fb"!. Me, = 300 
GeV, Te, = 2 GeV, Mq, = 50 GeV. All other effects 
have been switched off. 



ously: the MT2 "kink" method |19[, 
the invariant mass endpoint p6| , 27, 
28, 29 1 or the constrained kinematic 
method the polynomial intersec- 
tion method |2^], and p-r reconstruc- 
tion [pn, 31, 32], just to cite some pos- 



sibilities. Most methods, however ^ |2|, ||, ||], |2|l, require longer decay chains (at 
least two two-body decays) or a special topology, such as 4 on-shell intermediate resonances 
1 24] or 5 or more on-shell intermediate resonances [p^j. 

If each decay chain consists of a single two-body decay, where the visible one may or 
may not be a composite of visible particles, as described in Fig. |9|, one may use the MT2 
"kink method" ]19] to fix the value of mAr, i.e. exploit the fact that MT2™'^^ as a function 
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of the invisible particle trial mass has a "kink" for //jv = tun- The authors of |33| 
point out that in order to have a substantial change in the gradient ^^^^^^^^^—\vN=mN-: there 
must be substantial event-by-event changes, though. This can be triggered by substantial 
0{My) differences in the i^'tv, caused by the visible system being a collection of two or more 
particles, or by a large UTM. Otherwise the kinematics is so constrained that the gradients 
for /i^ < and > mj\i have to be the same, and no kink is possible. 

If a sizable UTM px is present, one may use the MT2_L method pi. This method 
uses the fact that N{iJ,]^), the number of times the MT2{ 

muN,pT) is larger than MT2{fii^,0), has a minimum for fj,^ = rriN. The advantage of 
using this method rather that MT2kink is that MT2(//jV) 0) may be calculated analytically 
and measured using the whole data sample, regardless the pt- This may be shown by 
using the fact that MT2{tiN,0) corresponds to MT2™^^ _L (^at), 

MT2??f'= min [max {M^j,_l, M^-p^}] , (5.8) 

F±=PlT±+P2T± 

the one-dimensional analogue of MT2, where 

M^^ll = m2 + + 2(^,p||^^|| - PiTii ■ pfT\\) 

and where _L and || refer to the projections of the px along the direction of the UTM. 
We will not discuss further the different methods to extract the two new particle masses. 



but refer the interested reader to the proceedings from the TeV 2009 |16| conference and 
to the review p3[| . We however wish to make a few remarks about the impact that valley 
radiation might have on these observables. Consider the MT2_L case, for example. 

In the presence of valley radiation one needs to consider two sources of deviations. 
Firstly, the tails coming from the interleaved radiation mechanism, see the MT2 distribu- 
tions in Fig. Secondly, as discussed in the previous section, in the presence of valley 
radiation, we expect the mean value of the invisible particle invariant mass to shift from 
its Breit-Wigner central value (M|^) = p.N towards the communicator mass Me^ (or Md^) 
value. We will show in subsection |6.2| that the MT2 distribution one obtains may be signif- 
icantly affected, see Fig. | for th LHC case. The MTl^f" (fiN) should be similarly affected. 
The number of events having MyifJ-Ni^T) > MyifJ-N^O) would then change accordingly, 
as would the minimum point fJ-N = ^N- 

We will return on the issue of the the trial mass and the radiation in subsection |6.2| . 
In the following, unless otherwise specified, the analysis will always assume = rriN. 



6. Effects of SUi?,)^, radiation at LHC 

At LHC the communicators are (mostly) pair-produced by gg or qq fusion and decay 
flavour diagonally into a SM d quark and a valley q^. For our study we assume the D^s to 
be spin 1/2 particles, and the q^s to be scalars. As earlier this choice affects the production 
cross section, but now both s- and t-channel exchange are involved, which complicates the 
pattern. Each radiates both SM gs and valley ^^.s. These in turn may radiate further gs 
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and QyS, respectively. Once the has decayed, the q radiates gluons, while the radiates 
g^s. The amount of hidden radiation emitted depends upon the valley coupling constant 
and on the mass ratio Mg^ /Md^ , see Fig. ^. At the LHC the communicator mass reach 
will be larger than for the ILC, so typically there will be more phase space available for the 
radiation. Both the valley gluons radiated by the and those radiated by the q^ have an 
impact on the visible particle distributions. The lighter the particle, the lower the cut-off 
scale for the radiation however, so it will be the qv that radiates the most, as before. 

When compared with the CLIC case, the LHC scenario presents several complications. 
Firstly one needs to convolute the production cross section with parton distribution func- 
tions. Thus the hard interaction — the production of the pair — no longer happens 
in or close to the center-of-mass rest frame. In this case it is crucial to consider longitudi- 
nally boost invariant observables such as MT2. Secondly, both initial- and final-state QCD 
radiation are more intense than the QED one is for the ILC case, resulting in a consider- 
ably larger upstream transverse momentum and an increased misassignment of radiation. 
Thirdly, there is an underlying-event activity that gives rise both to a generic low-p^ 
background and to occasional further hard partons that may be confused with the ones 
related to the valley process. Fourthly, the partons hadronize into more-or-less well-defined 
hadronic jets, the reconstruction of which introduces further smearing of the relevant kine- 
matic distributions. And finally, the set of possible background processes is much more 
varied and challenging to suppress. In our study we will take into account the first four 
points, but leave the last one to the experimental community, where already a large number 
of background-suppression techniques have been developed for various scenarios. 



6.1 LHC with 7 TeV 

The LHC will initially be running at of ^/s = 7 TeV and it is expected to deliver 1 fb^^ of 
data in 2010-2011. Under these conditions we need to consider much lower masses M^^ 
for the communicator than for the ultimate energy and luminosity case, in order to have 
large enough production cross sections, see Table |. 

Based on the above discussions, we choose to study the MT2 distribution, and specifi- 



cally its dependence on the v-radiation and on the a„ value. In Fig. |l^we have plotted the 
Oy dependence for three different mass values M^^ = 300, 400 and 500 GeV. The larger 
the mass difference M^^ — Mq^ the more phase space is available for the radiation. The 
smaller the Mq^, the lower the cut-off on the momenta, so the larger the amount of soft 
radiation. Given the low statistics, the (purely statistical) error bars on the endpoints 
are rather large, and even in the shoulder of the distribution it is hard to distinguish the 
curve with an = 0.1 from the curve with no radiation for the 300 GeV mass. In the 
intermediate case Md^ = 400 GeV, we need to have a rather strong ay = 0.2 coupling 
before the two curves can be separated. 



In Fig. [15| we show the MT2 dependence on the invisible q^ mass Mq^ , where the trial 
mass iiq^ is assumed to coincide with Mg^. Whether we look at the curves with a valley 
radiation (valley coupling ay = 0.1, left plot in Fig. |l^ or at the curves without hidden 



radiation (right plot in Fig. |T5|), the data points corresponding to Mq^ = 10,50 GeV are 
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Figure 14: The dependence of the MT2 distribution on the value for Md, = 300, 400, 500 GeV 
and Mq^_ = 10 GeV. The black curve corresponds to having no valley radiation. The red, green and 
blue curve correspond to a„ = 0.1, 0.2, 0.4 respectively. For the first 2 years LHC is assumed to run 
^/s = 7 TeV and to yield an integrated luminosity of 1 fb^^. The y axis corresponds to the number 
of events per 10 GeV mass bin, for this integrated luminosity of 1 fb"-'^. 



hardly distinguishable. The independence of the MT2 on the Mq^ value appears to be a 
characteristic which the radiation leaves unchanged. 
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Figure 15: Left: the Mq^ dependence of the MT2 when the valley radiation has a„ = 0.1. The 



blue curve corresponds to Mq^ =10 GeV and the red one to Mq^ 
dependence of the MT2 in the case of no valley radiation. Md^ = 
fb-1. 



= 50 GeV. Right: The same Mq^ 
400 GeV, y/s = 7 TeV and L = 1 
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Figure 16: Left: model a) decay — dqv with no v-radiation, Mq^ has a fixed value; right: 
model b) decay £)„ — )• dq^ with v-radiation, the system + g^s has an invariant mass distribution. 



One could argue that the fact that MT2 is hardly dependent on the Mq^ is due to the 
mass ratio M£)^/AIq^ being fairly large compared to the difference between the two Mq^ 
values we considered. Would this argument still hold true when the v-radiation is larger? 

For larger values of the one should consider the fact that the MT2 input parameter 
Mq„ corresponds to the mass of the g„ as seen by the visible particles. The visible particle 
momenta, i.e. the momenta of the d and the g radiated by it, which enter MT2, actually 
correspond to the invariant mass rather than the nominal mass value Mq^. 

In Fig. Q we show this invariant mass distribution and how this changes as a function 
of the coupling constant av As grows qv emits more and more valley gluons Qv The 
invariant mass of the + g^s system then grows, i.e the mean effective mass of the as 
seen from the SM q shifts from the bare Mq^ = 10 GeV towards M™*^^ = M^)^ . The energy 
which is left for the d quark then gets smaller as the (M|^) — )• M£,^ as shown by 

Ml. - 



where for simplicity we have put the = 0. 

Were this distribution a very narrow peak 
between the two cases in Fig. 16 so long as {Mqf^) = Mj'^, i.e. between a fixed M^^ value 



Were this distribution a very narrow peak around (Mqf^), there would be no difference 
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Figure 17: MT2 distribution function for the communicator masses Mo^ — 1 TeV when there is 
valley radiation a„ = 0.1 (red), 0.4 (green), 0.8 (blue) and when there is not (black curve), -y/s = 14 
TeV, Mq^ = 10 GeV. Number of events per 15 GeV bin. 



and an invariant mass distribution. One could speculate that so long as the mean value of 
the invariant mass (Mg^f) < Md^ the MT2 would basically remain unaffected. 

Imagine though having a v-radiation large enough to shift the substantially, 
M^^ — )■ M/)^, e.g. a = 0.4. One would naively think that replacing the Mq^ with (M^f) 
would give a better description of the case with radiation. However Fig. |^ shows that for 
> 0.1 the invariant mass distribution would also have a large spread around this central 
value. We will show in the next section that this spread causes further complications. It 
is precisely the spread in the distribution which constitutes the difference between case a) 
and case b), and it is this spread which is ultimately responsible for the different behaviour 
of the MT2 in the two cases. 



6.2 LHC with 14 TeV 

If we now assume that LHC will collect 100 fb^^ of data at center of mass energy ^/s =14 
TeV, then one may consider larger communicator masses, O(lTeV), and still deal with a 
sufficient number of events, see Table |2|. If the remains hght, the mass ratio M^^/Mg^ 
can be considerable and consequently the phase space available for radiation can be large. 
We expect the effects of the radiation to be significant. 



In Fig. 17 we show the dependence of the MT2 distribution on the hidden valley 
coupling constant , assuming we know the mass of the invisible particle so fig^ = Mq^ . 
In this case we are describing the strong coupling regime = 0.1,0.4,0.8. 

Notice how in this case one can actually separate (at least before any background or 
detector simulation is taken into account) the curve with = 0.1 and the one without the 
radiation. 

In the above study the assumption is that events are generated with a bare mass 
Mq^ and analyzed with the same (or almost the same) trial mass = Mq^. From the 
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Figure 18: Left: model a) comparison between different /i^^ = 1, 100, 350 GeV for Alq^ = 395 GeV. 
Right: model b) comparison between different ^iq^ for Af,„ = 10 GeV, a„ = 0.28, (M^f) = 395 
GeV. In both models M^i^, ~ 1 TeV, y/s — 14 TeV and the luminosity is assumed to be L = 100 
fb-i, 20 GeV bins. 



discussion in subsection 5.2, we could conclude that it is not very likely that we would 
know the mass of the Qy with high precision when valley radiation is present. 

In a hidden valley scenario the mass of the is assumed to be much lighter than the 
Dy mass, so we do not expect the MT2 to be very sensitive to mass differences of the 
order AMg^ < Mz),„. Indeed, the authors of |l|l and H show that the MT2™^^ mass 
dependence on the trial mass fiq^ is rather weak so long as < Mg^, whereas M.T2^^^ 
grows much more rapidly when ^g^ > Mg^. The fXg^ < Mg^ case is exactly what one 



observes in Fig. 18 



Let us however make the conservative assumption that we only know the order of 
magnitude of the Mg^ . Imagine trying to distinguish between the two models we described 
in Fig. a) the model with no radiation and b) the model with the radiation, when 
(Mqf^) = Mg^ . To be more concrete, assume a) has a fixed value mass Mg^ = 395 GeV 
and b) has an invisible particle mass Mj^^ = 10 GeV and = 0.28. We choose the value 
Oy = 0.28 so that {Mgf^) = M^^. 

In Fig. |l^, left side, we see the MT2 distribution for model a) for the case Md^ = 1 
TeV. As one may observe, this is a function of the trial mass /ig^, the best profile being 
the one with fig^ = Mg^. For ij.g^ <^ ^qv^ even for substantial changes in fig^ the MT2™'" 
does not change much. This essentially confirms what was reported by |jl^, 

On the right side of Fig. |l^ one may see the same distributions for the model with 
radiation. Contrary to what one would expect from the naive arguments given in the 
previous subsection, we see that choosing ~ (Mgf^) does not give the best description 
of the system. The MT2 curve overshoots the Md^ value by a good 10%. As anticipated 
in the previous subsection, this is due to the invariant mass distribution spread. Looking 
at Fig. 19, one may see that the invariant mass distribution has a wide spread, so event 
by event there could be large variations in the Mgf^. If one chooses a fig^ ^ (Mgf^) to 
analyze the set of events, for example the //^^ = 1 GeV chosen in Figure |l^, most of the 
events will have a "real" q^ mass, the invariant mass Mgf^, larger than the trial mass ^g^. 
Since MT2"^^^(/XgJ < MT2"^^(Mgf^) when < {Mgf^), these points do not contribute 
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Figure 19: The invariant mass distribution 
for the Qy and the trial masses we considered, 
= 1,100,350 GeV. 
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Figure 20: MT2 distribution when both qy 
invariant masses M^^°^ < fiq^ (upper) and 
when both are M^f^ > fiq^ (bottom), when 
= 10 GeV, ay = 0.28 and = 350 

GeV. 



to increase the MT2™'^^ value much, and the distribution will resemble rather closely the 
one obtains for = 10 GeV, apart from the softening in the shoulder. This is precisely 



what happens in the fiq^ = 1, 100 GeV curves in Fig. 18. 

When one takes a fig^ ~ {Mqf ) instead, e.g. fig^ = 350 GeV (which according to the 
naive arguments of the last subsection should have been the best of the three fiq^ choices), 
more and more events have a "real" invisible particle mass Mqf < fiq^. In Fig. |l^ these 
are the points to the left of the fig^ = 350 GeV line. These events will, if there is enough 
statistics, give MT2™''^(;UqJ >MT2™^^(Mgf '^), the real Md„ value. This is what happens 
to the curve for /ig^ = 350 GeV. 

To prove this point, we separately plot the MT2 distributions for the events with 
Mgf^ < fiq^ and the ones with Mgf^ > jjg^. As Fig. |2^ shows, for the former set 
MT2'^'^=^ > Md,. 

This said, we may return to the issue of distinguishing between the two models a) and 



b). As one may see in Figure 21, the two curves corresponding to the two models have 
very different shapes, the one with no radiation being the sharper one. Notice however 
that now the curve with radiation lies above the one without. This is not in contradiction 
with what we have shown in the previous sections. 

Summarizing, in the not-so-strong coupling regime ~ 0.1, when the invariant mass 
distribution is strongly peaked in Mq^, one may^ distinguish two models with the same 
values for Md^ and Mq^, one with v-radiation and one without. The curve corresponding 
to the no hidden radiation model will have a steeper drop in the shoulder region. 

For larger couplings, > 0.2, one should note that two MT2 distributions with the 
same endpoints may correspond to different M^^ values in the two cases, depending on 
whether the majority of the events have an invariant mass which is greater or smaller than 
the trial mass. A very conservative approach to solving this problem could simply be to 
assume = GeV or in any case fiq^ ^ (Mg^). 



^At least before background and detector analysis. 
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7. Conclusions 



We have here addressed the is- 
sue of detecting and identifying 
hidden radiation through its influ- 
ence on SM parton showers, and 
in general through its impact on 
visible particle kinematic distri- 
butions. We have done so in the 
context of Hidden Valley models, 
which we find are well suited to 
display the effects, but the specific 
models studied should be viewed 
only as representatives of a broad 
range of possible models with new 
symmetries. Thus, while we focus 
on the phenomenology of a fairly 
generic toy model, we also provide 
tools in Pythia 8 to simulate the 
effects of hidden radiation in various other hidden valley scenarios, e.g. different gauge 
groups, particle contents, and gauge and decay couplings. The novel feature in these tools 
is the interleaved SM and valley parton shower, i.e. the competition between visible and 
hidden radiation. 

Our preliminary study of the phenomenology of the toy model at e^e~ and at LHC 
colliders shows the following. 

At an 800 GeV ILC collider we could expect to observe hidden radiation for valley 
gauge couplings as small as > 0.05, so long as the mass of the communicator is smaller 
than 300 GeV and Mg^ < Me,- 

For the LHC phenomenology we need to distinguish between the first two year running 
with ^/s = 7 TeV and L = 1 fb~^ and later years with full design energy ^/s = 14 TeV and 
L = 100 fb-^ Whether one could observe hidden radiation or not depends strongly on the 
communicator mass, which determines the amount of statistics. The main signal of the 
hidden radiation — at least in our studies — is the softening of the shoulder of the MT2 
distribution. In the lower-energy case, for communicator masses around 500 GeV or higher, 
the statistics is so poor that one cannot expect to distinguish even a very strong coupling 
Oy = 0.4. For two models with and without hidden radiation, with equal M^^ masses 
in the [300, 400] GeV range and equal Mg^ < 50 GeV, one would need a valley coupling 
of the order a„ > 0.2 or larger to induce large enough effects on the MT2 distribution 
to distinguish between the two. In the higher-energy case, for order TeV communicator 
masses and Mq^ smaller than 100 GeV, the MT2 distributions show sizable changes already 
for = 0.1. 

We have also studied how the MT2 distribution depends upon the Mq^, the invariant 
mass and the trial mass In the case of a Hidden Valley scenario, Mq^ is always 
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Figure 21: Comparison between model a) Mq^ = 395 GeV 

and b) A/g\_ = 10 GeV, = 0.28, {Mj^i"^) = 395 GeV for 
= 350 GeV. In both models Md^ = 1 TeV, ^/i 14 
TeV and the luminosity is assumed to be i = 100 fb~^, 
20 GeV bins. Notice that in this case the curve with the 
radiation is the lower one. 
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assumed to be Mq^ <^ ■ Taking a trial mass fiq^, ^ M^^^ as an input parameter for the 
MT2 is thus a natural choice. When the mass is larger though, e.g. Mq^ ~ ^Dy, the 
issue of the trial mass Hq^ is no longer so trivial. The M^^ can be a broad distribution when 
radiation is present. Especially in the strong interaction case, ay > 0.2, {M^^) is strongly 
shifted towards M^^. This means that when one chooses a trial mass Hq^ ~ {M^^), roughly 
half of the events will have fiq^ > M^^ , causing the MT2™'^'^ to overshoot the real value. 
In this case the new masses thus have to be extracted from a combined fit, in which both 
masses and couplings enter as unknowns. 

Further studies of the background and detector simulations should follow, both for 
the kind of scenarios we have explored here and for other possible ones. This preliminary 
study, however, shows unequivocably that parton showers are a key tool in determining 
the presence of new hidden gauge groups and in the exploration of the hidden sector gauge 
group dynamics. 
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